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Abstract: The possibility that Lorentz symmetry is violated in gravitational processes 
is relatively unconstrained by experiment, in stark contrast with the level of accuracy to 
which Lorentz symmetry has been confirmed in the matter sector. One model of Lorentz 
violation in the gravitational sector is Einstein-aether theory, in which Lorentz symmetry 
is broken by giving a vacuum expectation value to a dynamical vector field. In this paper 
we analyse the effective theory for quantised gravitational and aether perturbations. We 
show that this theory possesses a controlled effective expansion within dimensional regu- 
larisation, that is, for any process there are a finite number of Feynman diagrams which 
will contribute to a given order of accuracy. We find that there is no log-running of the 
two-derivative phenomenological parameters, justifying the use of experimental constraints 
for these parameters obtained over many orders of magnitude in energy scale. Given the 
stringent experimental bounds on two-derivative Lorentz-violating operators, we estimate 
the size of matter Lorentz-violation which arises due to loop effects. This amounts to an 
estimation of the natural size of coefficients for Lorentz-violating dimension-six matter op- 
erators, which in turn can be used to obtain a new bound on the two-derivative parameters 
of this theory. 



1. Introduction 



Lorentz symmetry is widely considered to be one of the cornerstones of modern physics, con- 
firmed to an exquisite level of accuracy by experiment. These investigations are facilitated 
theoretically by the Standard Model Extension ]lj ^, which parametrises the renor- 
malisable Lorentz violating operators which can be introduced into the standard model 
Lagrangian. In this framework , some Lorentz violating perturbations are excluded at the 
accuracy of 10~ 27 . 

In light of this experimental evidence it is difficult to imagine that Lorentz symmetry 
could be violated at experimentally accessible energy scales. Additionally, theoretical dif- 
ficulties arise even if one considers energy scales which are presently inaccessible. Typical 
motivations for doing so include possible quantum gravity effects, for example, a Lorentz 
violating granularity of spacetime near the Planck scale. Recent arguments [0, ^ H sug- 
gest that effects of Planck scale physics can be communicated to low energy processes via 
radiative corrections, in particular, Lorentz violating contributions to low dimensional op- 
erators are argued to arise from the high- momentum part of loop integrals. Contrary to 
naive expectations based on the vast separation of scales involved, these contributions are 
not suppressed by inverse powers of the Planck mass and instead contribute at the percent 
level. This creates a new fine-tuning problem for these models. 

An altogether distinct possibility is that Lorentz symmetry is spontaneously broken at 



low energies, but is restricted to the gravitational sector [10, 11|. In this class of models, 
Lorentz symmetry is spontaneously broken by a vector field which is given a vacuum expec- 
tation value (VEV) using either a potential or a Lagrange multiplier constraint, examples 
of which include |J, g]| [lj, g]| ||, 0» © HI ® H ■ 



In this paper we consider a diffeomorphism invariant model, Einstein-aether theory [16 



17]. In this model, Lorentz symmetry is spontaneously broken by a dynamical vector 
field known as the aether field, A, with a timelike VEV, v, enforced through a Lagrange 
multiplier constraint, A. The Einstein-aether Lagrangian, in the mostly minus signature 
convention (+,—,—,—), and in units where h = c = 1, 

S = j ' d A x^[c 2 + C m + \[A 2 -v 2 )}, (1.1) 

where C m represents the matter content and 

C* = ~R- V^^V^V^, (1.2) 

where R is the Ricci scalar, K^ 1 ^^ 31 * 4 is the most general covariant polynomial of A 
consistent 2 with the Lagrange multiplier constraint, 

^-M1A»2M3M4 _ c ~^MlM2gM3M4 _|_ c ~gMlM3gM2M4 _|_ c ~ 3 gA»lM4gM2M3 _|_ fiiy^Ml gV3H4 ^ (1.3) 



1 See pi ||, ^] for a comprehensive review of experimental constraints. 

2 By consistent, we mean that there are no terms in which the number of A fields can be reduced 
covariantly by directly imposing the constraint. For example, any term proportional to A^VaA^ is not 
constraint consistent as can be seen by differentiating the constraint equation A 2 — v 2 . 
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and Cj with i = 1, . . . 4, are four free, dimensionless parameters. These are related to the 
usual Cj parameters, as found in p2| for example, by 



$ = (iVi» 2 Ci . (1.4) 

We omit the term R^A^A" because it is redundant; it is equal to a combination of the 
above terms under a field redefinition, up to a total derivative. Taking c% ~ 1, this model 
has two scales, the Lorentz violating aether VEV, v, and M which is related to the reduced 
Planck mass by, 

M 2 = il4 + ^±iV (1.5) 



as can be found in the weak field, slow- motion limit of the classical equations of motion [23]. 



In this paper we will work in the parameter regime v < M. For a review of the constraints 



on the c~j parameters see [22] and references therein. For some recent constraints arising 



from considerations of stability, see [24]. 

Implicit in the use of the two-derivative Lagrangian ( |1.2| ) in the solar-system or on 
larger scales is the assumption that the higher-derivative terms or quantum corrections are 
negligible. In this paper we make this statement precise using a dimensionally regulated 
quantum effective theory for Einstein-aether, making use of power-counting techniques. 

Using this formalism we are then able to investigate whether loop corrections can 
modify the scaling of the c% parameters. Recall that in QCD, for example, logarithmic 
corrections cause the value of the strong coupling to change from a s ~ 0.4 at 1 GeV to 
a s ~ 0.1 at 100 GeV, an effect which is subsequently measured p3]. If there is a scale 
dependence of the Cj in Einstein-aether theory, then it is of particular importance since 
the phenomenological couplings are compared with experiments which span many orders 



of magnitude in energy scale. For example, in the investigation of Cerenkov radiation [26], 
the most energetic cosmic ray particles considered have energies 10 11 GeV, which are many 
orders of magnitude greater than the typical curvature scales in the Solar-system. In turn, 
the length scales associated with the Solar-system are orders of magnitude smaller than 
those involved in cosmology. Nevertheless, observations on each of these scales have been 
used to constrain the two-derivative phenomenological parameters. 

Perhaps one of the most remarkable features of Einstein- Aether theory is the possibility 
that v ~ M, even in the face of otherwise stringent constraints. This large value for v is 
allowed because there are special combinations of the parameters, Cj, for which Einstein- 
aether theory is indistinguishable from general relativity in the PPN formalism The 
conditions on these parameters are independent of the size of v/M, allowing v ~ M. In 
this regime one has a large gravitational modification. If there is a log-running of the two- 
derivative parameters then it would be surprising if these parameter choices were preserved 
under a change in energy scale. In this case a value of (v/M) 2 < 10~ 7 would guarantee 
consistency with the PPN measurements without requiring delicately balanced parameters. 

This paper is organised as follows. We begin in section [2] with an introduction to 
quantum effective field theory techniques for gravitational theories, by reviewing the theory 
of quantised metric perturbations in pure gravity. We then apply this methodology to 
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Einstein-aether theory in section || where we define the effective Lagrangian for Einstein- 
aether, and discuss its perturbative quantisation. In section || we use power-counting 
techniques to estimate the size of a general Feynman diagram constructed in this theory, 
and in section [| we establish the hierarchy of a finite number of contributions to any given 
physical process. 

Having established that the effective theory possesses a controlled loop and derivative 
expansion we are then in a position to study the evolution of the parameters under the 
renormalisation group. We find that the low-energy, two-derivative parameters of this effec- 
tive theory do not receive logarithmic corrections. This result follows in a straightforward 
manner from the property that all interaction vertices are derivative couplings. 

Note that no steps in the above analysis could have been skipped; the non-running 
of the low-energy parameters follows from the structure of the perturbative low energy 
Lagrangian constructed in section ||, which in turn can only be trusted given the power- 
counting result (sections || and (5j) which demonstrates that there is a controlled expansion. 

In section ^ we discuss the size of the direct aether-matter couplings, and the associated 
potential fine-tuning problem. Assuming a sufficiently small or vanishing two-derivative 
direct coupling term, as measured, we estimate the expected size of Lorentz violation which 
arises due to loop effects. This amounts to an estimation of the natural size for coupling 
constants of four-derivative dimension-six operators, for which there are also strong exper- 
imental constraints. This size depends on the ratio of scales v/M, hence it can be used to 
constrain the two-derivative Q-parameters. We also comment on the use of a momentum- 
space cut-off regulator. 

Related works on the subject include the semi-classical approaches [28, 29] which begin 
with a direct coupling of quantum matter to classical aether and metric, and then see 
how loops of matter renormalise the aether theory. Our interest is different here, we 
wish to study effects associated with quantised metric and aether perturbations. The 
phenomenology of tree-level diagrams in Einstein-aether theory has also been considered 
l2l 



2. General relativity as an effective field theory 

In this section we provide a brief introduction to effective field theory by considering the 
quantum theory of metric perturbations in pure-gravity ]30| , 31]. The Einstein-Hilbert 
Lagrangian contains a heavy mass scale, the (reduced) Planck mass, 

Co = -^-V=gR. (2.1) 

Perturbing this Lagrangian about a weakly curved background solution g^ u = g^ u + 
Mp 1 ^^, it is clear that canonically normalised metric perturbations carry inverse powers 
of the heavy mass. The graviton operator of lowest dimension is dimension four, the kinetic 
term, whilst all other operators are interacting and therefore carry increasing powers of 
1/Mp. These are irrelevant operators, signalling that the theory is non-renormalisable. 
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The presence of only non-renormalisable interactions means that divergences arising 
from loop integrals cannot be absorbed by redefining the couplings present in the La- 
grangian. However, provided we are only interested in low energy scales, E <C Mp, then 
effective field theory offers a consistent solution to this problem. Using a regularisation 
procedure which preserves the symmetries of the Lagrangian, loop divergences are guaran- 
teed to have a counterterm if one adds all irrelevant local operators consistent with these 
symmetries. 

Adding all irrelevant operators consistent with diffeomorphism invariance gives an 
effective action with an infinite number of coupling constants, which we have arranged in 
a derivative expansion 3 

S = j d A x^j (-^R + axR^R^ + a 2 R 2 +..X (2.2) 

The meaning of renormalisability in the effective field theory sense can then be seen most 
explicitly in the background field method, where explicit calculation confirms that one- 



loop |32|] and two- loop [33] divergences are of the form of these higher derivative terms. 



Indeed, the effective theory for gravity is fully renormalisable in this sense [34]. 

Of course there are now an infinite number of terms and an infinite number of unknown 
coupling constants. Therefore if this theory is to be of any use one must find a consistent 
way to truncate the Lagrangian. To do so one must make a further assumption - that 
the coupling constants are natural. Here this means that all coefficients are 0(1) up to 
combinatorial factors and the scale Mp. Of course the size of the cosmological constant 
term is at odds with this assumption - a problem which is not resolved by the use of effective 
field theory. With this assumption, one can perform a straightforward power counting for 
the size of a general diagram Ai computed using a dimensional regulator and constructed 
with t loops, E^ external graviton legs and iVw vertices with d derivatives, 



E X 2+2£+E d=2 A r(d) (°!-2) 



M ~ Mp- E -< (i- ) , (2.3) 



where E is the typical energy scale constructed invariantly from the external momentum in 
the graph, E <C Mp. The hierarchy of contributions to this process may then be read off 



from (2.3). The leading contribution to a particular process are tree level diagrams with 
two derivatives. Next-to-leading contributions are one-loop diagrams constructed with two 
derivative operators and tree level diagrams containing one four-derivative operator, and 
so on. The expansion ( |2.2[ ) is therefore an expansion in (E/Mp) 2 , which breaks down if E 
approaches Mp. Now we can see the utility of the effective approach; whilst we have an 
infinite number of undetermined coupling constants, this does not imply that the theory 
is not predictive. To compute a process of typical energy E, one simply has to specify 
the desired accuracy of the result, which determines where one can truncate the derivative 
expansion. This leaves a finite number of parameters, rendering the theory predictive. 



3 Note that we have not included the super-renormalisable cosmological constant term. This is determined 
by experiment to be incredibly small ~ y/— g 10~ 120 Mp, and can be neglected if one is not interested in 
cosmological scales p5|. 
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Note that one is only able to draw the above conclusions because there is one scale in 
the problem, Mp. This is a challenge presented by Einstein- aether, where the physics of 
the two-derivative Lagrangian depends on two scales. 

In summary, the effective field theory of perturbative quantum gravity has the advan- 
tage of allowing quantum corrections to be calculated, provided the energy scales involved 
satisfy E <C Mp. At these energy scales, the theory is renormalisable as an effective field 
theory @. 



3. Einstein-aether effective theory 

In this section we construct an effective theory of quantised metric and aether perturbations 
for Einstein-aether. Specifically, we use a path-integral formulation with the Lagrange 
multiplier constraint directly imposed. This construction allows for the computation of 
quantum effects and their subsequent renormalisation. 

We begin by writing down the most general action consistent with diffeomorphism 
invariance and the Lagrange multiplier constraint. It will be convenient to use a dimen- 
sionless aether field in what follows, A = A/v, where v is the size of the VEV for A. 
Throughout this paper any field with a tilde is canonically normalised. To represent the 
matter content of this theory we have added a scalar, (ft, for simplicity. We have, 

S = J d i x^{C[A^g tlv ,<ft\ + X{A 2 -l)), (3.1) 
C[A», gil¥ ,4>] = E «£a, s )°5 a ' S) K'^^], (3.2) 

d,a,s,i 

where Q^j denotes a constraint consistent, diffeomorphism invariant term which can be 
constructed from d derivatives, a aether fields, s matter fields and the metric. There are 
potentially a number of different terms with the same counting due to differing Lorentz 
structures, and so we add a further label, i, which enumerates all possible terms. Note 
that we have not excluded the potentially problematic direct (ft, A^ couplings, for example, 
A^d^(ftA u d u (ft; we have included all terms. 

As in the pure-gravity example, there are an infinite number of terms in the Lagrangian 
( |3.2| ). We therefore need a way to truncate this theory, and this will be the focus of sections 
^| and ||. In this section we will define the perturbations and the quantisation procedure 
using a dimensional regulator. The main reason that we use a dimensional regulator is 



to preserve the diffeomorphism symmetry. This regulator is used in pure-gravity [32, 33] 



and in theories which are structurally similar to Einstein-aether in the absence of gravity, 



for example in chiral perturbation theory [35]. Dimensional regularisation allows for a 



controlled truncation of the effective theory as we shall show, whilst a momentum-space 
cut-off gives an uncontrolled effective expansion. We discuss this issue in section |(| 
The coupling constants have mass dimension, 



(i) 

(d,a,s) 



4-d-s, (3.3) 
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and with the assumption of naturalness are to be built from the scales in this theory, 
together with some combinatorial factors and an 0(1) number. The scales that appear in 
these couplings are dependent on the renormalisation point \i at which the theory is defined. 
Here for convenience we will focus on the characteristic energy scale associated with the 
curvatures in the solar-system, or on larger length scales. Hence it is most convenient to 
consider the effective theory defined at a renormalisation point well below the mass of the 
lightest particle in the standard model, which we denote by m. We take, 

£~/i<m<t><M. (3.4) 

That is, we have 'integrated out' all massive particles through a matching procedure 4 on to 
the effective Lagrangian fl3.2| ). The results obtained within this regime are easily extended 
to higher energy scales, m < E <C v < M , and we discuss this in section ||. 

At low energies only virtual processes of the massive particles are important and this 
is exactly what the higher-derivative operators in the effective Lagrangian (|3.2| ) are chosen 
to achieve through the procedure of matching. Hence they will typically have inverse 
powers of the scale m in their coefficients. Also important at low energies are of course 



the massless degrees of freedom. The two-derivative Lagrangian (1.2) explicitly defines the 
massless gravitational and aether modes. Additionally we have the massless remains of the 
standard model. For simplicity we will idealise this matter content by using the scalar field 
above, 4>, taking it to be massless and to not self-interact at two-derivative order. 

We perform a 3 + 1 split of A, and a decomposition of the metric about flat space, 

A = ad + TT i d i , (3.5) 
g fJiU dx fJ/ dx u = (rjnv + j^dx^dx" , (3.6) 

Treating the ir and 7 fields perturbatively picks out one effective field theory in particular, 
corresponding to choosing a coordinate frame in which the background value of the aether 
field is purely time aligned. In principle one may choose to treat a different set of fields as 
perturbations, corresponding to perturbation theory about a different background Lorentz 
frame, and a different effective field theory 5 . However the frame we have chosen here is the 
one in which the phenomenology of Einstein-aether has been investigated and compared 
to experiment, for example in the PPN formalism [38, 27]. 

The approach we take here is to eliminate a from the Lagrangian using the Lagrange 
multiplier constraint, which may be performed formally in the path-integral by integrating 
over A, 

Z = j VAV 1 V\ e iJd 4 xC+\(A 2 -i) x j VAVl6 (42 _ ^ e ifd*xC_ ( 3J ) 



4 For a discussion of continuum effective field theory see for example [ [36[ |37[ 

5 It has been recently pointed out |24| that if the observer frame is highly boosted with respect to the 
aether rest frame, higher derivative terms can be lower order in the energy expansion than naively expected. 
Again, we stress that this is a different effective field theory altogether, and one would have to perform a 
new power counting analysis. 
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This path-integral requires gauge fixing, and there is the same amount of diffeomorphism 
degeneracy here as in the case of general relativity. One may therefore use a pure- 
gravitational gauge breaking term of choice, such as the Prentki gauge 6 . A purely gravi- 
tational gauge fixing of course leads to the introduction of the usual Faddeev- Popov ghost 
terms, and consequently, graviton-ghost interaction vertices. Being quadratic in the ghost 
field, these ghost terms have the same power counting properties as pure graviton vertices. 
Hence, as is the case for power counting in pure gravity (see for example [01]]), we do not 
need to include them explicitly in our analysis. One should therefore bear in mind that 
in this paper we have omitted Feynman diagrams which contain internal ghost lines and 
these would have to be included if one went beyond power counting. 

Continuing from (|3.7| ) we may integrate over one of the fields in the delta-functional. 
In particular we choose to integrate over A = a. Employing the identity 

J Va Det (^p) & {A 2 - l) = 1, (3.8) 

and the expression for a given by solving the constraint equation A 2 = 1, 

a = %Q (-90iir l ± [500 + (goiir 1 ) 2 - 5oo5ij7rV] , (3.9) 

we obtain, 

Z oc J VttV 7 [Det (±2( 5o0 + (goi^f - 5oo^7rV)i)] ~* e'/^^WI, (3.10) 

where the plus and minus roots correspond to a choice of future or past-directed A, respec- 
tively. There will of course be Goldstones associated with the vacuum Lorentz violation, 
and these correspond to a mixture of the new tt fields and the metric perturbations. There 
are no massive modes, since the constraint forces the aether to lie exactly on the unit 
hyperboloid. 

The determinant may be absorbed into the Lagrangian by the use of the functional 
identity DetM = expTrlogM. This will result in new, quartically divergent operators 
in the n fields. An analogous determinant appears in the perturbative quantisation of the 
non-linear sigma model, see for example [[K], 41 1. In this case the determinant is known 
to account for the non-linearly realised symmetry, retaining the correct invariance of the 
functional integration measure. In turn, it gives rise to new power-law divergent operators 
in the Lagrangian which cancel against loop integrals which do not respect the symmetry 
of the Lagrangian. Here we work with dimensional regularisation for which there are no 
power-law divergences, and the determinant will not contribute new interaction terms. If 
a cutoff regulator were to be used instead, we expect that the determinant here plays an 
analogous role to the determinant found in the non-linear sigma model, and such terms 
must therefore be included. 

All operators in this theory are the result of expanding covariant, constraint consistent 
terms. As an example, consider a kinetic term for the aether field in the 2-derivative La- 
grangian. Ignoring the index structure and order-one numbers, its perturbative expansion 



'See for example the lectures by Veltman in 
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will be of the form, 

v 7 — 5 cv 2 VAVA ~ c ydftdir + (jj^ dydir + {^j^j djdj + -TtdiidiT + -^irdirdy + . . , 

(3.11) 

where the perturbations have been canonically normalised, 

7T = VTT, 7 EE Mj. (3.12) 

Thus once the integration over A and a has been performed, one covariant term in the 
Lagrangian gives rise to a finite number of marginal operators and an infinite tower of ir- 
relevant operators of increasing mass dimension, with one coupling constant up to factors of 
v/M. Then, if our regularisation and renormalisation procedure preserves the symmetries 
of the covariant form of the Lagrangian, then all logarithmic divergences may be absorbed 
by these couplings. For example, the effective theory of pure gravity discussed in section || 
has been shown to be renormalisable in this sense @. We do not have an analogous proof 
for Einstein-aether, and we assume this property. 

Note that in general the propagator for the 7 and ir perturbations is not block- 
diagonalised, that is, 

+ 0. (3.13) 

This of course depends on the specific numerical relationships between the coefficients of 
the marginal operators, and the chosen gauge fixing condition. For the purposes of power 
counting it will be convenient to retain the perturbations as defined here, because we are 
interested in counting the mass scales v and M. 



4. Power Counting 

As discussed in section an effective theory with an infinite number of coupling constants - 
such as that outlined in the last section - must be supported by a power-counting argument 
which allows for a controlled truncation of the Lagrangian based on the accuracy of the 
result required. This is necessary for the theory to be predictive. In this section we take the 
first step; we perform power-counting to estimate the size of any Feynman diagram. This 
is supplemented by section [|, where we show that this leads to a truncatable Lagrangian. 

We begin by considering a general Feynman diagram in momentum space with t loops 
and En,E^ and E^ external aether, graviton and matter lines respectively. We denote its 
amplitude by A4, which has the mass dimension 

[M] = 4 — E w — E 7 — E$. (4.1) 

As usual we will not include factors of the momentum arising from external leg propagators 
or the momentum-conservation delta function in M . Of course, if one wanted the size of the 
(Ej,+ E^ +i£y)-point function itself, then these would have to be included in addition to the 
amplitude M. Since we are working in dimensional regularisation any renormalisation scale 
fi will not contribute to the dimension of Ai. The dimension ( ^4.1[) is therefore accounted 
for entirely by typical light external momenta E, and the heavy mass scales v, M and m. 
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Calculating M by combining contributions from vertices and internal lines generically 
results in a complicated multi-integral expression. However an estimate for its size is easily 
deduced since the leading powers of v,M,m necessarily factor out. Our first step will be 
to count these heavy mass factors, which are all contributed by couplings at vertices and 
by the non-canonical internal propagators. 

As previously stressed, the sizes of the higher derivative couplings, and some of the 
two-derivative couplings are not known. Ultimately these are to be determined by ex- 
periment. We first discuss the two-derivative couplings, for which the theory is currently 
over-parametrised, since we may always perform a simultaneous rescaling of the aether field 
and redefinition of v to absorb one of the couplings. There are four possible terms which 
can yield a kinetic term for the spatial aether perturbation, it. These are, 



(4.2) 



a Um v^V^V^A,, a (2)4 , ) V^A^A v V„A p V v AP. 
We use the ir kinetic term to define the scale v, that is, we take 

"(2,2,0) ~ "(2,4,0) ~ ( 4 - 3 ) 

Additionally the Ricci scalar term defines the scale M, and we canonically normalise <f>. 

Hence for the two-derivative Lagrangian the only terms not yet discussed are direct 
couplings between the aether and matter. The matter Lagrangian here is massless and 
does not self- interact and so is invariant under a constant shift, <j) — > <fr + c. Demanding 
that this shift symmetry is not broken by direct aether matter couplings then restricts the 
derivatives to act on the cj) fields and not on the aether fields in these terms. Hence the 
only two-derivative direct aether-matter coupling with these symmetries contains exactly 
two aether fields. This term contains a Lorentz-violating contribution to the <j> propagator 
at leading order in its perturbative expansion, 

"(2,2,2) V ZI gA> i A v d^d^ = a (2)2i2) d Q 0d o + • • • (4.4) 

where the ellipses denote higher dimension operators. Given the experimental bounds on 
the smallness of direct aether- matter coupling constants (see for example 0, ||, ||], ||, |(| ) , we 
impose the minimal phenomenological requirement that this marginal term is not larger 
than the contribution from the Lorentz-invariant piece of the <ft propagator. That is, we 
require 

"(2,2,2) < 1- (4-5) 

This concludes a discussion of the size of the two-derivative Lagrangian couplings. 

Moving on to the higher-derivative couplings, we include all operators which are con- 
sistent with the required symmetries, using the dimensionless fields g^u, A and the unit 
mass-dimension cf). Our goal here is to establish whether or not the theory has a controlled 
effective expansion, and to achieve this we construct these couplings with inverse powers 
of the mass scale m, 

a (d,a,s) ~ m(d - 4)+8 ' ^ 4 > rn<v<M. (4.6) 
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We choose these values because they constitute an upper limit on their natural sizes, since 
any loop corrections they receive will be smaller by factors of m/v or m/M. If these 
couplings are indeed smaller than (|4.6|), as may be found through experiment, this will 
only improve the effective expansion. 

We now look at how mass scales enter in the internal propagators. In momentum 
space, the aether and metric marginal operators in the Lagrangian take the following form, 

I lip) \ T I M 2 A flu v 2 B^ u \ ( 7(-p) 



PfiPu, 



(4.7) 



where 7 is a vector of the ten components of the metric perturbations, 7? is a vector 
of the three spatial aether components and T denotes the transpose on these indices. 
The matrices A M!y ,B^ and C^ u account for order-one numbers and the structure of the 
marginal operators in terms of the ten components of 7 and the three components of tt. 
The spacetime indices of these matrices contain the structure of the momentum dependence 
of these operators, and it is worth emphasising that they are not Lorentz tensors. To 
determine the typical sizes for the inverse of this matrix we take the limit where v/M is 
small. Then, 



M 2 A^ v 2 B^ u 



u 



A"' 

M B 



U' 



Ml 



^(C- 1 B T A- 1 )..,. + 0(^ 



M 



1 

w 



IF 



v_\2 
M) 



u 



M 



1 
M 1 



1 

W 2 



C-WA- 



where, 



U 



o 



(A-^C- 1 ) +0 



Mi 



TP 



. M 2 



Ml 
vl 



(4.8) 



(4.9) 



(C _1 B T A _1> 



In general, the off-diagonal blocks of this matrix, (A _1 BC _1 ) and _ . , 
will not vanish. Then, the heavy mass scales entering into the propagator sizes for the 
non-canonically normalised perturbations tt and 7 will be given by, 



1 



7' 



7 



TT ■ 



TT ■ 



■ TT 



7 



M 2 p 2 
1 

v 2p2 



1 



M 2 p 2 



1 + 0| M 



(4.10) 



In the special case where v ~ M, there will be only one scale, and hence the above 
propagators also give the correct counting in this regime, despite their derivation in the 
regime where v <C M. Hence we assume that they provide the correct counting for v < M. 
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With a different choice of gauge for the perturbations it is possible that the off-diagonal 
blocks could be made to vanish. However, any physical results will be independent of the 
gauge fixing condition, and so we are free to work in the generic case for which the above 
counting holds. Here the p 2 factor represents the contribution of the typical light scale 
of the momentum in the propagator, and in general this will not be a Lorentz invariant 
quantity. Note however that the momentum dependence is still quadratic 7 in both the 
time, and the spatial components of p. 

With the mass scales contributed by the propagators determined in ( 4.1 Q|) , one could 



explicitly count the mass-scale contribution from the (ivy) propagator simply by counting 
how many times it appears in the graph. One must be careful however, since when tt 
and 7 are both canonically normalised, the propagator (fry) still carries a factor of v/M. 
Ordinarily, specifying the external leg content determines completely which propagator 
is used in the external leg. However, here we must account for the possibility that the 
'species' can change in the external leg, which subsequently contributes further factors of 
the ratio v/M. A useful mnemonic to keep track of the factors of v and M introduced 
by the {ttj) propagator is to count it as a 'vertex' which has one 7 leg and one tt leg, 
and which carries a the coupling constant v 2 p 2 . We emphasise that this is not a Feynman 
vertex, but a combinatorial tool. Hence it is subject to the additional rule that it may be 
included only once in any given line of the Feynman diagram. We then have, 

1 1 

V p* M z jr 

as required. 

We enumerate the total number of vertices by N, which we then partition according 
to the number of derivatives, d, the number of matter fields, s, and the number of factors 
of the aether field, a, present in the original covariant term from which the operator is 
derived, 

N = ^ AtW<m) + v> ( 412 ) 

d,a,s, 

where V counts the number of species-changing 'vertices', as described above. We denote 
the number of internal propagators in our diagram by I^,I n and 1^. 

In the general case, summing up contributions from heavy mass scales we obtain, 

{ "> (4 13) 

x m ^>^^ a ' S) ^-ds) E D v -E^ M -E^ 

The remaining mass dimension of M is accounted for by D powers of the light scale E, 
and the last two terms account for the non-canonical normalisation of 7 and tt. We may 
now determine D by the requirement that Ai have the correct mass dimension (|4.1| ) . After 
some manipulation, and use of the topological identity, i~ + I n + 1^ = £ + N — 1 we find, 

D = 2 + 2£+J2 N(d ' a ' S) (d-2)- (4.14) 



7 This is to be contrasted with other theories which manifestly break Lorentz invariance such as the ghost 
condensate model f42| where dispersion relations become quartic. 
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There are now no unknown quantities in ( 4.13 ), and we have therefore estimated the size 
of any Feynman diagram in this effective theory. 

The equation ( [4.14 ) is the usual effective field theory result found by Weinberg in 
the context of Chiral Lagrangians [43|. It is the same as the ^-counting result in the 
pure gravity case ( |2 . 3[ ) . Increasing the number of loops in the diagram or the number of 
derivatives at vertices necessarily increases the power of E, which, in a theory with only 
one mass scale has the simple interpretation of increasing suppression by factors of E/M. 
However, such a result is meaningless in the presence of extra scales if one cannot show 
that there is a controlled expansion in ratios of those scales too, for example, v/M. We 
will address this issue in the next section. 

As an example of the counting performed in this section, consider the following dia- 
gram which we have constructed using the interaction vertices: y/— gX7AX7A D y 2 7r 2 e?7T(?7r, 

V 2 "f"fdjdTT, 




(4.15) 

This graph has five vertices in total, N = 5, which are made up from two interaction 
vertices from the two-derivative Lagrangian A^ 2 ' 2,0 ) = 2 (filled circles), and three species 
changing 'vertices', V = 3 (crosses). The internal line content is I& = 0, I % = 2, J 7 = 3 



and the external line content is E^ = E^ = 0, E n = 4. Inserting these values into ( 4.13] ) 



we estimate the size of this diagram to be, 



M 



-) (-)' 
M \MJ 



(4.16) 



5. Validity of the effective expansion and logarithmic corrections 



Using the power-counting results of the last section we show that the effective theory 
defined in section || can be truncated to a Lagrangian with a finite number of terms, with 
the number of terms dependent on the order of accuracy required. As discussed, this is 
necessary for the theory to be predictive. Once this is established we are then able to 
estimate the scale-dependence of the couplings arising due to logarithmic corrections. 

In order to determine whether there is a controlled expansion we consider the size 
of a general diagram ( |4.13j ). One can perform a straightforward re-arrangement of the 
summation ( |4.13| ) to obtain, 



M ~ m 



' E~ 


D 


"772" 


m 




. V . 






' V ' 




X 








lm\ 



2U+E 7r +2I^+E^-2N < - 2 - 2 '°1 -27V( 2 A°) -2Ar< 2 > <°) -2V 



27 7+J B 7 -2iV( 2 ' . ) ^(2,2,2) 
[«(2,2,2)J 



(5.1) 



All quantities in square brackets in this expression are 0(1) or smaller, and hence there is 
a well defined expansion if the power to which they are raised is bounded from below. If 
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they are not bounded from below then they can provide large ratios of scales which could 
offset the E/m suppression arising from increasing loop or derivative order. In this case 
diagrams from arbitrarily high derivative order can contribute and the derivative expansion 
cannot be truncated. 

It is easily seen that the power of «(2,2,2) i s bounded from below. For powers of v/M 
we note that the total number of graviton 'ends' in the graph (2I 7 + E^) cannot be less 
than the total graviton valency contributed by the Ricci-scalar term, which is greater than 
two per vertex, that is, 

(5.2) 



2/ 7 + £ 7 > 2iV (2 ' ' 0) . 



The remaining power of m/v is also bounded from below, as can be seen by a similar argu- 
ment; each vertex in the graph arising from the expansion of \J— gVAVA , \J— gAASJ A\J A 
or y/—gR, must have a combined total graviton and tt valency of at least two. Therefore 
the total graviton and tt valency in the graph cannot be less than this amount, that is, 



24 + E w + 2 J 7 + E 1 > 2 (iV(2,2,0) + at(2A0) + #(2,0,0) + i/j _ 



(5.3) 



We have now established the existence of a controlled expansion. That is, given the 
energy scale of interest and the order of accuracy required, there are a finite number Feyn- 
man diagrams to evaluate. This analysis justifies the use of the two-derivative Lagrangian 
in describing processes at sufficiently low energies. We remind the reader that in order 
to obtain the power-counting result ([[0]), we constructed the higher derivative couplings 
using inverse powers of the scale m. The coupling constants as determined experimentally 
could be smaller, in which case only vertex contributions are affected and the effective 
expansion would be improved. 



We are now in a position to discuss the full implications of equation ( 4.14 ), in the 
context of the renormalisation group. A consequence of is that a one-loop graph 

necessarily scales as at least E 4 . For example, consider the matrix element for tt — tt 
scattering up to order E 4 , which we denote by 



M(E) = Y J M i 



(5.4) 



where i runs over all graphs which can contribute. We work in a subtraction scheme, 
for example MS. Once any divergences from loop integrals have been absorbed into the 
couplings, suppressing 0(1) numbers and any Lorentz structure, we find 



M(E) 



E 2 



+E 4 



Yj b j a (2,2i,0)0) 



3=1 



4 , 
7 = 1 



(5.5) 



+ (E 6 ) 



where bj and b y *' are dimensionless 0(1) constants, and the sum runs over the number 
of aether fields in the covariant terms which can contribute to this process. The maximum 
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number of aether-fields which can occur in a covariant term containing d-derivatives is 2d, 
as any more cannot be constraint consistent. A is a constant of mass dimension [A] = —4, 
which depends on details of the size of the higher-derivative couplings. For the assignment 
of higher-derivative couplings used above, A ~ v ~ 2 m~ 2 , but in practise it can take a smaller 
value and the theory will still have a good expansion. The expression for M{E) must be 
independent of the renormalisation point, fi, 

^-M(E) = 0. (5.6) 

Hence in general the four derivative couplings evolve under the renormalisation group. 
The two derivative couplings do not receive radiative corrections 8 , and hence we conclude 
that the two-derivative parameters do not have any scale dependence. This justifies the 
application of experimental constraints to these parameters, despite their spanning of many 
orders of magnitude in energy scale. 

We now comment on the regime in which the above results were obtained. One can 
work in the effective theory above the mass scale m, for which there is at least one massive 
particle which can contribute to the running of the two-derivative gravitational parameters. 
It is easily seen that the only possible one-loop corrections to the two-derivative gravita- 
tional/aether parameters involve at least one mass term and hence are at least (m/M) 2 
suppressed, instead of (E/M) 2 as in the massless case. These corrections are therefore 
negligible for a standard- model matter content. 

Of course there are very heavy masses in the solar-system, the Sun being one example. 
The quantum effective theory for non-relativistic masses in the case of General Relativity 



has been studied in 44}] . Despite appearances, there is still a good effective expansion 
for this theory, which is more clearly seen when one adopts a non-covariant, physical 
gauge. The expansion parameter was found to be the same as in the vacuum case, that is, 
(E/Mp) 2 . It is unclear whether Einstein-aether theory possesses a controlled expansion in 
the presence of non-relativistic masses and we leave this as an open question. 

6. Direct aether-matter couplings 

In this section we use the effective theory to estimate the expected size of Lorentz violation 
due to loop effects. This is intimately related to the size of the coefficients of direct- 
aether matter couplings, and we outline a way one can use these estimates to bound the 
two-derivative q parameters. 

The two-derivative direct aether-matter couplings are commonly set to zero (for the 
matter content above, this coupling is a(2,2,2))- This is due to strict constraints coming 
from the observed accuracy of Lorentz invariance in the matter sector, together with the 



The order E piece of equation (5.5) together with the condition ([jjj), suggests that it is possible for 
the two-derivative coupling constants to evolve on some hyper-surface in parameter space. This would mean 
that the two-derivative Lagrangian was over parametrised. However, the Lorentz structures associated with 



the different contributions Q(2,2,o) an d «(2,4,o) are m general different, and this is not captured by (5.5). 
Furthermore, if this invariance of the Lagrangian exists then it would be seen classically, and no such 
symmetry of the two-derivative Lagrangian has been found. 
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fact that this coupling modifies the matter dispersion relation in a Lorentz violating way 
(see [PI ). 

Here, to make progress, we assume that there is some property of the underlying 
theory - a symmetry perhaps - which predicts a sufficiently small or zero direct coupling 
at two-derivatives in this low energy theory. Supersymmetry provides one such candidate 
P§| ) |46|| . In section [5] we saw that the parameters of the two-derivative Lagrangian do not 
receive radiative corrections when below the scale m and only small corrections above, and 
so the approximate size of this coupling does not change with the renormalisation point 
in this theory. Hence by setting this coupling to zero at one renormalisation point we will 
find it remains zero at larger length scales. Of course there can still be Lorentz-violating 
effects in matter processes which arise through loop effects, and higher derivative terms. 
We also know that the higher derivative terms can be scale dependent and so in general it 
is inconsistent to set these to zero. 

Hence we allow for higher derivative direct couplings, for which there are significant 
experimental constraints in the four-derivative case, though these are less severe than for the 
two-derivative operators. For example, recent work for fermionic operators [47] constrain 
these kinds of dimension-six couplings to be, a < W~ e M~ 2 . For simplicity we retain the 
use of scalar matter, <j>. At four derivative order we have two kinds of direct-coupling term 
which, once expanded, give rise to a dimension-six operator for eft alone. These are, 



cn42p)y/=gAAVd<l>Vd<l>, a^^y/^gAAAAVd^Vdet), (6.1) 

where we have suppressed Lorentz indices. 

To estimate the size of these couplings we use naive dimensional analysis (NDA). This 
method was first proposed in |^|] for the chiral theory of low energy QCD, which provides 
a way to estimate geometrical factors of 4ir in the non-renormalisable couplings. In this 
section we use the NDA approach to estimate factors of the scales v and M. 

The basic idea behind NDA is that the coupling constants are set to be the same size 
as the largest radiative correction they each receive. The NDA result then gives the typical 
size of the coupling which one would generate under a change in the renormalisation point, 
H. 

To study radiative corrections to 0(4,2,2) an d 0(4,4,2) we are f ree to study any process 
to which they contribute at tree-level. Since we are using power-counting to determine 
the size of these diagrams, any cancellations which occur between individual loop graphs 
cannot be seen. Power counting therefore gives the best estimate for the corrections to 
these couplings when we study the process for which the leading loop correction is the 
smallest. If larger corrections are found by studying other processes, covariance of the 
unperturbed Lagrangian then guarantees that there will be appropriate cancellations. 

The processes which give the best estimate for the size of loop corrections to (^O]) are 
those which involve at least one external tt leg. Here we consider the process E n = = 2, 
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E~ = 0, for which we can construct the following Feynman graphs, 




E 4 

a (4,a,2)-^T, 




E 4 ( v\1 
a {i,a,2)^T {m) ! 




E 4 ( v \^ 



(6.2) 



Hence, the contributions to this process coming from vertices obtained in the perturbative 
expansion of the covariant terms ( |6.lD are governed by the leading J 7 = diagrams. 

Radiative corrections to these coupling constants come from one-particle irreducible 
loop diagrams with the same number of external legs and the same power-law E scaling 
(we saw in section || that loop corrections with a different E scaling correct coefficients of 
vertices containing a different number of derivatives). Hence, by (4.14) these can only be 
one-loop diagrams constructed from the two-derivative Lagrangian. Using (|5.1|) we find 
that these diagrams have the size, 



M 



loop 



E\ /m\ 4 f v \ 
m J 



2U, 



v) KMJ 



(6.3) 



Finding the largest loop correction therefore amounts to minimising the power of v/M 
in ( |6.3| ). It is easily seen that without a direct aether- matter coupling constant at two- 
derivative order, a one-loop graph with both external tt and <ft legs must contain at least two 
internal gravitons, J 7 > 2. Possible loop graphs with J 7 < 2 are either not connected and 
therefore do not contribute, or not one-particle irreducible. Hence the leading contribution 
is given by I 7 = 2. Examples of these graphs include, 




(6.4) 



Hence to order E we have the following contributions to the E^ 
process, 



E, 



2, E~ 



M(E) 



E q 



a 



(4,a,2)(A t ) + 



1 



v \- 



KM) l ° g 



E 



+ 0(E 6 ). 



(6.5) 



Again, having subtracted any divergences (in MS, for example) and then using the property 
(|5.6j ), we see that under a change in the renormalisation point an order correction to 
a (4,a,2) (/^) 1S generated. Hence, as an NDA estimate of the size of these couplings we obtain, 

«(4,a,2) 



V 

M 



M 2 ' 



a 



2,4. 



(6.6) 



These couplings are therefore naturally smaller than the usual (E/M) 2 suppression by a 
factor of (v/M) 2 . Hence if v is sufficiently small compared to M, then M ~ Mp and 
the dimension-six Lorentz violating <f> operators are naturally small. An exact calculation 
of these corrections to a more realistic matter content could be used as a constraint on 
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the size of v/M, given the experimental bounds found in [47], for example. This would 
subsequently provide an upper limit on the size of the usual phenomenological parameters, 
Cj = (v/M) 2 Ci ~ (v/M) 2 . Tentatively assuming a similar bound for our idealised matter 
content we obtain 9 , a < 10~ 6 . 

We will now comment on the use of a momentum-space cut-off regulator. Even if one 
were to overlook the requirement of symmetry preservation in the renormalisation proce- 
dure, there would still be the problem that this cut-off regulated theory has no controlled 
expansion. This is because there are potentially an infinite number of one-loop corrections 
from vertices containing arbitrarily many derivatives, all of which could contribute the 
same order of correction. In this set-up, in contrast with the dimensionally regulated case, 
the couplings in the two derivative Lagrangian could receive power-law corrections, the size 
of each individual correction dependent on the value of the cut-off chosen. For example, 
taking A ~ v and using the NDA approach to estimate the size of the two-derivative direct 
coupling, one finds 



i: 



■l 



"(2,2,2) ~ {jj) • (6-7) 

Whilst it is interesting that this suggests the theory does not require finely-tuned cancel- 
lations if v/M is small enough, it is inconclusive given that there is no controlled effective 
expansion. 

7. Discussion 

In this paper we considered the effective theory for quantised aether and metric pertur- 
bations of Einstein-aether. With the inclusion of scalar matter we estimated the size of a 
general Feynman diagram, and demonstrated that the theory can be truncated in a con- 
trolled way when a dimensional regulator is used. We showed that the phenomenological 
Q parameters of Einstein-aether theory receive only negligible logarithmic corrections, jus- 
tifying the scale independence which is assumed when they are compared to experiment. 
We estimated the expected size of Lorentz violation for matter due to loop effects commu- 
nicating the Lorentz violation between the aether and matter sectors, and outlined how 
this could be used as a new constraint on the Cj parameters. 

With a controlled way to truncate this effective theory it becomes possible for more 
detailed features to be investigated. It would be interesting to study the inclusion of 
non-relativistic masses into this low-energy theory, and to see whether a good effective 



expansion is retained. This has been investigated in the case of general relativity in [44], 
where a good expansion is found. Furthermore, the leading quantum corrections to the 
Newtonian potential has been calculated for pure gravity p9| , |30| . It would be interesting 
to see how the aether perturbations modify this result, and its dependence on the observer 
frame in which the perturbations were defined. 



9 There is a potential caveat here. If one invokes supersymmetry at higher energies in order to suppress 
any relevant and marginal Lorentz violating operators, then any E 4 contributions to the matter dispersion 
relation are naturally small ^] . In this case the coupling constants for higher- dimension supersymmetric 
operators are naturally compatible with experimental constraints, and hence a weaker bound would be 
obtained. 
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There are of course a number of technical issues which have not been addressed here. 
A step towards a background field method for this theory may be found in the context of 
chiral perturbation theory [^] , which provides a way to deal with a similar type of length- 
fixing constraint in a covariant way. In pure gravity infrared divergences can be cancelled 
by graphs with additional soft gravitons [50], something which has not been investigated 
here. Additionally, perturbation theory in pure gravity is renormalisable in the effective 
field theory sense |34j|; for Einstein-aether we have assumed this property. 
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